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We study the superfluid to Mott insulator transition of a mixture of heavy bosons and light 
fermions loaded in an optical lattice. We focus on the effect of the light fermions on the dynamics 
of the heavy bosons. It is shown that, when the lattice potential is sufficiently deep to confine 
the bosons to one dimension but allowing the fermions to freely move in three dimensions (i.e. a 
mixed-dimensionality lattice), the fermions act as an ohmic bath for bosons leading to screening and 
dissipation effects on the bosons. Using a perturbative renormalization-group analysis, it is shown 
that the fermion-induced dissipative effects have no appreciable impact on the transition from the 
superffuid to the Mott-insulator state at integer filling. On the other hand, dissipative effects are 
found to be very important in the half-filled case near the critical point. In this case, in the presence 
of a finite incommensurability that destabilizes the Mott phase, the bosons can still be localized by 
virtue of dissipative effects. 



I. INTRODUCTION 

The interest in systems of interacting bosons and 
fermions has been a recurrent and central topic in the 
study of the many-body problem. Many early studies 
were concerned with dilute solutions of "^He in ''He (see 
e.g. [l], for a review) as well as with the problem of elec- 




FIG. 1: Schematic representation of the system studied in this 
work: A mixture of light fermions and heavy bosons loaded 
in an anisotropic optical lattice. As the bosons are assumed 
to be heavier, the are confined to one dimension by the lattice 
potential. 



trons coupled to phonons in solids (see e.g. [2]). This 
research led to the understanding of important phenom- 
ena such like the polaron and Cooper pairing [2 . More 
recently, these concepts have reemerged in the context 
of ultracold atomic gases H] , where new types of in- 
teracting Bose-Fermi mixtures have been experimentally 
realized [51ll6|. Indeed, such experiments with ultra- 
cold gases have made it possible to study and envision 
Bose-Fermi systems P"7H5S] that can exhibit very differ- 
ent properties from their condensed-matter counterparts. 

Thus far, much research has focused on understanding 
how interactions with the bosonic component of the mix- 
ture influences the properties of the fermions and, in par- 
ticular, how the interactions mediated by the bosons can 
possibly induce fermion superfluidity (see e.g. [31 [311 [31] 
and references therein). The complementary problem, 
namely, understanding how the properties of bosons are 
modified by their interaction with fermions in a mixture 
has only recently attracted interest, especially motivated 
by a series of ground-breaking experiments with Bose- 
Fermi mixtures loaded in optical lattices [BlfTTl ITS] . 

Within this setup, in recent years a number of groups 
have addressed the problem of how the addition of 
fermions to a Bose gas in an optical lattice affects the 
phase transition from superfluid to Mott insulator in the 
latter 7, 9 11,15, 24,25,^- Thus, experimental obser- 
vations have been reported indicating that fermions ef- 
fectively decrease the quantum coherence of the bosons. 
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thus making it easier for the latter to become Mott insu- 
lating [3 [9HTT| [15] . In the case of attractive interactions 
between bosons and fermions, this effect has been ex- 
plained by a 'self-trapping' effect: the bosons move in the 
lattice potential that effectively becomes deeper by the 
addition of fermions to which the bosons are attracted. 
This 'self-trapping' would have the opposite effect in the 
case of repulsive boson-fermion interactions. However, in 
this case, the two components were found not mix in a 
deep lattice [TT] . 

In connection with the experiments referred to above, 
there has been some theoretical discussion on other, per- 
haps more subtle, effects of adding fermions to an in- 
teracting boson system [24j |25]. These effects concern 
the physics of polarons, where one particle (in this case, 
the bosons) is dressed by its interactions with a different 
species (the fermions). Thus, the boson-boson interac- 
tions are screened, becoming much less repulsive. The 
bosons also undergo dissipative effects, which involve the 
creation of real particle-hole pairs and other kinds of low- 
energy excitations in the Fermi gas. Indeed, as discussed 
below, within a weak coupling approach the 'self-trapping 
effect' arises at first order in the strength of the boson- 
fermion interactions, whereas the polaronic and dissipa- 
tive effects arise at second order. 

Going beyond mean-field theory, Yang [ID] has stud- 
ied the effect of the boson-fermions interaction on the 
superfluid (SF) to Mott insulator (MI) transition in a 
three dimensional Bose-Fermi mixture. He found that 
the properties of the transition at the particle-hole sym- 
metric point (i.e. at the tip of the 'Mott lobes') is mod- 
ified and becomes either a first order transition or a dif- 
ferent (i.e. not XY) second order transition. However, 
in this work we find that in a mixed dimensionality sys- 
tem, the universality class of the transition (2D XY) is 
not modified by the boson-fermion coupling for integer 
filling. The latter only introduces screening of the exter- 
nal periodic potential and the boson-boson interactions. 
Whereas the former tends to make the system less (more) 
superfluid for attractive (repulsive) boson-fermion inter- 
actions, the screening of the boson interactions always 
favors superfluidity. On the other hand, for a half-filled 
lattice, we find that the quantum phase transition from 
the SF to the charge-density wave (CDW) phase is modi- 
fied by the presence of the fermions. However, transition 
remains continuous and belongs to the 2D XY universal- 
ity class. 

The outline of this article is as follows. In the following 
section, we introduce the basic model of a Bose-Fermi 
mixture that will be subsequently analyzed. There we 
also outline the derivation of its effective low-energy de- 
scription. In Sect. |III| we consider the effect of the Fermi 
gas on the Mott insulator to superfluid transition of a 
Bose gas confined to one dimension [HJ ITT^HS] . The per- 
turbative renormalization group is used to analyze the 
low-energy properties of the effective low-energy model, 
the effect of the fermions on the superfluid to insulator 
quantum phase transition is studied. Depending on the 



lattice filling and the boson-boson interactions, the Mott 
insulator can be stabilized at integer or half-integer fill- 
ing, and the effects of the Fermi gas are very different on 
both transitions. Thus, we have separated the discussion 
into two subsections, III A and III B| Finally, in Sect. [V] 
the main conclusions of this work are summarized. 



II. BASIC MODEL 

A. Hamiltonian 

The system under study is an ultracold mixture of 
bosonic atoms (mass ms) and single-species fermionic 
atoms (mass mp) loaded in an optical lattice (see Fig. [I]). 
The repulsive interaction between bosons is described by 
an interaction potential vbb{t^ — r')- The latter can be 
either the Lee- Yang-Huang pseudo-potential, which ac- 
counts for the s-wave scattering of ultracold atoms (as 
for alkali or alkaline earth atoms), or a dipolar potential 
(as for Chromium or polar bosonic molecules). Further- 
more, fermions and bosons are assumed to interact only 
via a short-range potential, which is also described by the 
Lee- Yang-Huang pseudo-potential. Inter-fermion inter- 
actions are negligible because, by the Pauli principle, the 
dominant scattering channel for single-species fermions 
is p-wave, which, away from resonances, is very weak 
at ultracold temperatures. The optical lattice potential 
UBiF){r) = f/J^r^ sin2 kLX + U^j^'^sm^ fc^y + sin^ fc^z), 
where — 2t:/Xl and Al is the laser wavelength. 
It is further assumed that ^ Uq^, that is, the 
bosons move in a strongly anisotropic two-dimensional 
lattice. We further assume that the bosons are heavier 
(i.e. uip/rriB ^ 1), which means that their motion along 
two directions (here y and z) is strongly suppressed be- 
yond zero-point motion, thus effectively confining them 
to one dimension for at least the duration of the exper- 
iment. However, the fermions, being lighter, can hop 
in all three dimensions but the large laser intensity re- 
quired to create the strong confining lattice potential for 
the bosons, which implies that Uqj^ <C U^^^, the fermion 
dispersion will be anisotropic (see below Eq. |6]). Thus, 
the Hamiltonian reads: 



H — Hb + Hp + Hbp 

r 



H 



B 



dr 



2mB 



V4'|j(r)V*s(r) + t/B(r)/5B(r) 



I 



+ - / dv' pBir)VBB{r-r)pBir') 



Hp 



dr 



-2mF 

Hbf = 9bf I dv pb{v)pp{v), 



V*Hi-)V*F(r) + C/F(r)/5F(r) 



(1) 
(2) 

(3) 
(4) 
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where ^b(F){^) is the boson (fermion) field op- 



'B{F) 

erator, which obeys 



yr),#s(r' 
- r)) 



(5(r - r) 

(anti-)commuting oth- 
the boson 



erwise; PBiF)ir) = *b(f) W 
(fermion) density operator and A^b(f) = / dr pb(f){'^) 
the boson (fermion) number operator. The boson- 
fermion interaction is parametrized by the couphng 
gBF = 2Trh?aBF /Mbf, where Mbf — mBmp j [tjib + 
mp) is the reduced mass and a^F is the s-wave scatter- 
ing length. Since we are interested in the ground state 
phase diagram in the thermodynamic limit of the above 
system, we have neglected the harmonic trapping poten- 
tial, which is also present in the experiments. Note that 
an implicit assumption of our analysis below is that the 
bosons and fermions are mixed. For short range interac- 
tions between the bosons (i.e. for Vbb{'^) — 9bb5{v)) the 
problem of the bosons and fermions forming a uniform 
mixture in the lattice geometry studied here has been 
previously considered in Ref . [32] . One conclusion of this 
work is that the uniform mixed phase in this Bose-Fermi 
system is always stable provided the density of bosons 
and fermions is sufficiently high, for both attractive and 
repulsive interactions (see Ref. ^3? for further details). 

The Hamiltonian introduced in equations ([I] [2j |3| , and 
Q contains too much information about energy scales 
in which we are not interested. Since our goal is to an- 
alyze the ground state and low-lying excitations of the 
system, we next derive an effective Hamiltonian that is 
much more appropriate to this end. The first step is 
to project the Bose and Fermi fields onto the lowest 
Bloch band of the lattice potential. Thus, we expand 
^-B(r) - Sr^^o(i"J- -R) '^B-R.{x) where wo(r_L -R) are 
the Wannier orbitals describing particles localized round 
the site R = ^{m,n)\L of a 2D (square) lattice. For 

the fermions, *I'_F(r) — X^k ¥'k(r) /k, where </5k(r) are the 
Bloch states of the lowest band. Note the differences in 
treatment of the Bose and Fermi fields, which refiects 
their differences in mobility introduced by the conditions 
discussed above. Hence, upon neglecting terms coupling 
different lattice sites, the bosons are described by 









j dx 


2mB 


da^-^Bnix) 





R 

f i ^ y dxdx' Vbb{x - x')pbr{x)pbr{x'). 

Ft 



However, the fermions are described by: 



(5) 



(6) 



where the sum is over k belonging to the first 



BriouUin zone and e(k) = e||(/s) + ej^(kj^) _ 2„j. 

2tj_ (cos kybo + cos fcz^o): where bo = is the lattice pa- 
rameter, and we have assumed that the periodic potential 



along the x direction is so weak that efi^ectively amounts 
to a renormalization of the fermion mass. Finally, the 
boson-fermion interactions are described by: 



Hbf — gBF 



R 



dr \wo{rj_ - R)|^ pBR{x)pF{r), (7) 



where r = {x,y,z) = (a;,rj^). In the above expres- 
sionm we have approximated the boson density operator 
PB{-r) = Pb{x, r_L) Er |w^o(r_L - R)|^ Pbii{x). 



B. Integrating out the fermions 

The total Hamiltonian obtained upon projection onto 
the lowest Bloch band H = Hb + Hp + Hbf is still 
too complicated to solve. Since we are mainly interested 
on the low-temperature properties of the heavier bosons, 
which are much slower, a first step towards understanding 
the latter is integrating out the fermion degrees of free- 
dom. To this end, we rely on the path integral representa- 
tion of the partition function Z = Tr e-''[^-'"«^«~^^^^l 
for the Hamiltonian, H = Hb + Hp + HpB , which allows 
us to write: 



Z = [dipBd^'BdipFdipp] e' 



(8) 



where 

S = Sb 
Sb 



Sp + Sbf, 

^ J dx J dTlljBR{x,T)dr1pBR{x,T) 



R 

Pb 
h 



^ J dx J dT\lpBR{x,T)\'^ 

Hb{t). 



R 

dr 

^ 

hp 



(9) 



Sf = J2 L drfix, k) drfik, t) - ^/(k, r) 







P-F 



Sbf — 



1 

dTHBp{T). 



(10) 

(11) 



where (3 — {kBT)~^ is the inverse of absolute tempera- 
ture and kB is Boltzmann's constant. Thus, the effective 
action for the bosons is defined by the following equation: 



-Sotf [-i/J-B ,-0-E 



[dfdf] 



-S B—S F—S BF 



(12) 



where {. . .) p = Tt: pp . . . and pp = Zp^ p-fs{HF-t^Nj.) ^ 
being Zp = Tr e^^'-^^^^^^-' the non-interacting fermion 
partition function. To make further progress, we shall 
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assume that the interaction between the bosons and the 
fermions is pertubatively smah. Therefore, the above 
functional integral can be performed using the cumulant 
expansion, which yields: 

— p-{Sbf) + ^{S%p-{Sbf)'^) + --- 



(13) 



The leading term is 
{Sbf)f 



9BF 

h 



^ / dr I dr|wo(r_L 



R) 



X pBR(x,r)p^(r), 



(14) 



where p%{v) — {pf{y,t))f is the equilibrium density 
of the Fermi gas (in the absence of the bosons). Since 
p%'{v) is periodic, (14) amounts to a correction to the 
periodic potential that the boson gas undergoes. The 
correction has the same sign as the coupling qbf, which 
means that e.g. for attractive boson-fermion interactions, 
the effective potential seen by the bosons is deepened by 
its (mean-field) interaction with the fermions. This ef- 
fect has been termed 'self-trapping' and has been stud- 
ied both theoretically [15] and experimentally . We 
shall not study it any further here. Instead, we focus on 
the second order term, which leads to much more inter- 
esting physics. Neglecting the coupling between different 
sites R (i.e. terms where R' ^ R) yields: 

-\{SIf-{Sbf?) = ^^T.I drdrdr'dT'pBii{x,r) 



R 



xxf{x 



t')pbr{x',t'), 



(15) 

where r — (x,rj^), r' = (a;',r^). After defining 
Fo(r_L,r_L) = |u;o(r_L)i«o(r^)P, wc introduce 

Xf{x,t) = y (ir_L«ir^i^o(r_L,rl)XF(r,r',T), 



XF(r,r',r) 



:{5pF{r,T)SpF{r\0))i 



(16) 



Thus, up to 0{gBF): '^'^ obtain the following effective 
action for the bosons: 



cff.R 



hl3 



dr 



hp 







dr i dx 



dx i^*BRix, T)dri^BTi{x, t) 

h 



2m j 



-\dxipBR{x,T)\' 



+ 



1 



^^dT 

h 

h 



dx 



Ub\\{x)-p.b |V'SR(a;,r)| 



+ -y / dxdx'\il}B-R.{x,T)\'^ 



X Vbb{x - x')\'4>BTl{x' ,t)\'^ 

+ ^ J dxdTdx'dT'\^B-Rix,T)\'' 

X Xf{x, x', t - t')\'iPbti{x' , r')p. 



(17) 



where Ub\\{x) = UB\\{x)+gBF J dr± \wo{r±)\'^ p°p{x,r±). 
Note that we have thus reduced the problem to a set 
of one dimensional systems independently coupled to a 
fermionic bath. Therefore, in what follows we shall drop 
the lattice index R and study the phase diagram of a 
generic ID system coupled to the fermionic bath. 

However, one important caveat is in order when con- 



sidering the applicability of the effective action, Eq. (17 1. 
In what follows, we will not treat the bosons and the 
fermions on equal footing. Such a treatment would re- 
quire to also account for the effect of the bosons on the 
fermionic component of the mixtures, which may mod- 
ify the density response XF(r, r',r). Nevertheless, below 
we shall assume that xf(i", r', ''') is well described by the 
non-interacting limit where we take gsF = 0. Indeed, 
this assumption is qualitatively correct as long as the 
Fermi component of the mixture remains a Fermi liquid, 
which is reasonable given that the fermions are much 
lighter, interact with the bosons weakly, and therefore 
their energy is dominated by the kinetic energy. How- 
ever, strictly speaking the bosons will mediate effective 
fermion-fermion interactions, which, at sufficiently low 
temperature, lead to a pairing instability of the Fermi 
gas. Since the gas contains a single species of fermions, 
such a paring instability takes place in a high angular mo- 
mentum wave (most likely, p-wave) and at relatively low 
temperatures compared to the Fermi energy iip- Given 
that present cooling techniques in optical lattices cannot 
reach temperatures below a few percent of pF, we can 
safely neglect this possibility. Other instabilities that can 
gap the fermion spectrum, such as a charge density wave, 
occur at particular values of the lattice filling and/or lat- 
tice parameters and we will also neglect them in what 
follows. 



C. Low-energy effective theory 



In order to deal with the effective boson model in 



Eq. (17), we shall use the method of bosonization [HI 144) . 
Thus, we first integrate the high-energy density and 
phase fluctuations of the bosons, and introduce two col- 
lective fields, 6{x), and (l){x) describing the phase fluctu- 
ations in each ID system. In terms of these fields, the 
Bose field and density operators read: 



^Bix)^Apl/^ e^'^-\ 

Pb{x) = ^^Bi^)^B{x) 



Po 



-dx(t>{x) 



-I- Po ^ B'm cos 2m [(j){x) + kpx) 

m>0 



(18) 



(19) 



where po = Nb / (AIL) is the linear density of bosons 
in each of the M ID systems of length L of the lattice 
and kp — Trpo- The amplitudes A and Bm depend on 
the microscopic details of the model and can not be ob- 
tained using bosonization. Using the above expressions 
and retaining only the most relevant operators in the 
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renormalization-group sense yields 



The above action, Eq. (20), provides an effective de- 



1 



27tK 



dx 



- idr<j)f + V {dx 

V 



(20) 
(21) 



S'„U1 = - ^ / drdx cos {2p(t){x, t) + xSj,) , (22) 
7ra5 J 

where we have introduced the foUowing notation: v is the 
sound velocity of the ID Bose gas whereas vK = Vp ~ 
hkp/mB and K/v is proportional to the system com- 
pressibility; flo ~ hv/iiB is a short-distance cut-off j44j . 

The bare dimensionless coupling of the term describ- 
ing the periodic potential in the bosonization language 
is gu — UQ\\{BpPoa'^) I2hv . The cosine term with p — I 
describes the effect of the potential in the case of integer 
fiUing of the lattice, 6p^i = 2{kL — kp) being a mea- 
sure of the incommensurability of the system. However, 
near half filling, we must consider the p = 2 term with 
dp=2 = 2fci — Akp as a measure of the incommensurabil- 
ity. In the half-filled case, the above effective Hamilto- 
nian describes the transition from a Tomonaga-Luttinger 
liquid (TLL) to a fractional Mott insulating state which is 
also known as a charge density wave (CDW). The stabil- 
ity of the CDW state requires smaller values of the Lut- 
tinger parameter K than those that are achievable in 
the Lieb-Liniger model [55] describing bosons interacting 
via a short range potential in a ID wave-guide [S5], for 
which the minimum value of K is one |44| . Smaller val- 
ues of K are accessible when either the bosons posses a 
dipolar moment [Ml |50] or in the so-called super- Tonks 
regime [iil[5T1 [55] . 



FIG. 2: In the limit of strong one-dimensional (ID) con- 
finement, the inter-tube coupUng can be neglected and the 
system of Fig.[T]can be mapped to a single ID tube immersed 
in a Fermi gas. The Fermi gas has two important eflfects: it 
introduces screening of the boson-boson interactions and it 
behaves as a dissipative bath that introduces quantum dissi- 
pation. Assuming that a periodic potential is applied longi- 
tudinally to the tube(s), which drives a transition from the 
superfluid to a Mott insulating phase, we study the modifica- 
tion of the phase diagram due to these effects. 



scription of the low-temperature properties of the boson 
system which includes (through the renormalization of 
the potential Ub\\ Ub\\) the effect of the Fermi gas at 
the mean-field level. The dynamical effect of the fermions 
on the bosons is taken into account, to leading order in 
5bF: by the last term in Eq. (17). However, since the 



dynamics of the (heavier) bosons described by (21) is 
much slower than the lighter fermions, some further sim- 
plifications of (15) are possible. First, we note (see Ap- 
pendix|B]) that, at T = 0, the fermion density correlation 
function introduced above, xf{x, x' , r) can be written as 
follows: 



Xf{x - x',t) 



+ 00 



duj 



.\r\ 



lmx§{x-x',u), (23) 



where Xpi^ ~ 2;',w) is the retarded version of the same 
correlation function. We have also assumed, consistently 
with what was stated above, that the effect of the peri- 
odic potential can be neglected as far as the calculation 
of xf{x — x',uj) ~ Xpix — x' ,uj). The above expression 
allows us to treat separately the high frequency density 
fluctuations from the low-frequency fluctuations of the 
fermionic gas. This can be done by introducing the fol- 
lowing response functions: 

X<(x,r)=/ — 5(a;)e-Hlmxf(x,a;), (24) 



xUx,r)= I — 5c(w)e-"l"llmxf(x,c^), (25) 



where g{uj) is a frequency cut-off function, which can be 
chosen in various ways as the result will be largely in- 
dependent of this function; 5c(w) = 1 — g{(^)- Below we 
use g{uj) = e~'^'^<=, where <C max{^, ^}. The cut-off 

frequency ~ is chosen such that the high-frequency 
density fluctuations of the Fermi gas can adapt instan- 
taneously to the (slow) dynamics of the boson density 
fluctuations described by pb{x,t) (cf. Eq. 19). Thus, 



dxdrdx dr pb{x, t)xp{x — x ,t — t )pb{x , t ) 



= / dxdx'dtdr pb{x,t + ^)xpix ~ x',t)pBix,T ~ ^) 



I 



dxdx' dT pb{x, t)xp{x — x' ,ui = 0)pb{x' , r) 

= I dxdx' dr pb{x,t)xf{x — x' ^ Q)pb{x\t) 

dxdx'drpBix, t)xp{x — x' ,lu = Q)pb{x' , r), 



(26) 



(27) 

where x^(x — = 0) = / dtxp{x — x' ,t) and similar 
definitions for x'f{x — x' ,uj = 0) and xf{x — x' ,ui = 0). 
Therefore, the effective action describing the interactions 
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between the bosons mediated by the fermi gas takes the 
form: 



5, 



cS,BF 



9%F 



dxdx'dr pb{xtt)xf{x — x' ^ 0) 



2h 

y- Pb{x' ,t) + ^^^^ j dxdx'drdr' pb{x,t) 



(28) 



X r(a; — x'.T ~ t')pb{x' , r'), 

where the dissipative kernel T{x,x',t) is defined as: 

T{x - x', t) = r) -xU^- x\io - O)J(t). (29) 

Note that, by definition, J dr r{x—x', r) — 0. This kernel 
can be evaluated as follows. Since we assume the Fermi 
component of the mixture to be a Fermi liquid, we note 
that for the latter —lmxp{x — x', uj) (x u ioi w ^ \^p\ [T]. 
In the present system, the small lo limit of this function 
is obtained explicitly in Appendix [B] at T — 0. It can be 
written as 



lmxF{x — x' ,L0 ^ — ) = —TrD(x — x')u). 



(30) 



where D{x) is a positive function of x which is com- 
puted in Appendix [Ej Introducing this expression into 
(pll), yields: 



r(x — a;', t) 



D{x-x') 
{\r\+r,? 



(31) 



at T = 0. Introducing the above expression into Eq. (28), 
we arrive at: 



5, 



cS,BF 



Af 



dxdx'dr pb{x,t)xf{x — x' ,Ll! ~ 0) 



2h 

X pb{x' ,t) — ^^-^^ J dxdx'drdr' pb{x,t) 



D{x~x') 
\t~t'\+Tc) 



-Pb{x',t'), 



(32) 



The results of the model calculation described in Ap- 
pendix|B]for the functions D{q)/h = — Im x|?((7, w) (for 
w <C and the static response function Xpi'i^^ — 0) 
are displayed in Figs. |4] and [3| It can be seen that both 
functions are rather smooth (i.e. non-singular) functions 
of the longitudinal wavevector q. This assumption will 
prove important below. Furthermore, for certain values 
of the lattice filling, which determines the Fermi energy 
e^, see Appendix |Bj D{q) can be made negligible or zero 
for wide ranges of the wavevector q. This opens the pos- 
sibility of tuning the strength of the dissipative effects by 
simply changing the fermion density. Note, however, that 
by strongly reducing the fermion density, the stability of 
the mixture may be jeopardized 

Thus we see that the boson interaction mediated by 
the Fermi gas consists, at low frequencies, of an instanta- 
neous part (which stems for high frequency density fluc- 
tuations of the Fermi gas) and a dissipative part, which 
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FIG. 3: The static response function Xf(9i^ ~ 0)/^ of 
the fermions for different values of the Fermi energy ep = 
2.0, 1.5, 1.0, 0.5, -0.5 in units where = 1 and W = 4t_L = 
1. See Appendix [B] for details of the calculation. 
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takes the form of a retarded ^ \ interaction. The latter 
stems from the excitation by the motion of the bosons 
of real low-energy particle-hole pairs, which in a Fermi 
liquid yield the linear-o; behavior of the density response 
function (i.e. Landau damping). As discussed above, the 
instantaneous part of the interaction can related to the 
static density response of the Fermi gas and leads to a 
renormalization of the sound velocity v and Luttinger pa- 
rameter K describing the low-temperature properties of 
ID boson system. The renormalized parameters obey: 



v{gBF) v{g_ 



BF 



0) 



K{gBF) K{gBF = 0) 



2^XF(g = 0,c. = 0), 



(33) 



Furthermore, since the fermion-induced interaction is a 



density-density interaction (cf . first term in Eq. 32 ) , we 
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FIG. 4: Imaginary part of the fermion response function 
divided by the excitation frequency, cu, for — >■ 0'*', for dif- 
ferent values of the Fermi energy ep ~ 2.0, 1.5, 1.0, 0.5, —0.5. 
Units where — = 1 and W — 4ti = 1 have been used. See 
Appendix |B] for details of the calculation. 
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have that [H] : 

v{gBF)K{gBF) = v{gBF = 0)KigBF = 0). (34) 

These equations describe, to lowest order in gsF, the 
screening of the boson-boson interaction by the fermion 
gas, which leads to corrections to the parameters K and 
V in Eq. (21), which depend only on the boson-boson 
interaction. 



Using the bosonization formula (19), we obtain the 
representation of the dissipative action in terms of the 
density field 4>{x,t): 



— 

qb _ 



9fD 



dLxdrdr 



{\r-T'\+T,Y 



(35) 
(36) 



d.drdr' ^°^'['^("'"|''^(!'"')] . (37) 
ao J {\t-t'\+t,Y 



In the derivation of the above perturbations to the Gaus- 



sian action, Eq. (21), we have retained only terms whose 
integrands are not oscillatory and are the leading terms 
in a gradient expansion. However, in the case of a half- 
filled lattice, the following term: 

Sl = -'-^ I dxdrdr' y "I + '^^f; , (38) 

must be also taken into account. This dissipative umk- 
lapp interaction arises from the periodicity of the boson 
system, for which which at half- filling Akp = is a 
reciprocal lattice wave number. In this regard, we must 
recall that, in a periodic system, the (lattice) momentum 
along the x direction is conserved modulo a reciprocal 
lattice wave number. Note that this term will be also 
generated by the reno rma lization gro up f low from prod- 

(cf. IstI. 



uct of the Su (cf. Eq. I22j) and 5^ 



Furthermore, the bare dimensionless couplings are: 



9d{0) 


= 5|f^(? = 0). 


(39) 


9bDiO) 


= 2glpByaoD{q = 2A:|), 


(40) 


9ud{Q) 


= 2glpByaaDiq = 2kf). 


(41) 



In the above expressions we have made explicit the de- 
pendence of the couplings of the cut-off scale, oq through 
the parameter £ = log that is ao{i) = e^ap and thus 
£ = corresponds to the scale of the bare cut-off ag ~ vTc, 
being Tc the short-time cut-off introduced earlier. 



III. RENORMALIZATION GROUP ANALYSIS 

Physically, the renormalization group (RG) flow of a 
system describes its behavior as it is cooled down to- 
wards the absolute zero. The effect of temperature can 
be mimicked by decreasing the short wavelength cutoff 
introduced to properly define the low-temperature 



effective model of the last section. As the absolute tem- 
perature decreases, the ground state is approached, and 
the couplings that define the effective low-energy theory 
of equations (21 1, (22), (36) (37), etc. (i.e., K, v, gu, gbD, 
...), must change accordingly in order to account for the 
reduction of the available excited states. Thus, the quan- 
tum phases of the system can be studied by analyzing the 
asymptotic behavior of the 'fiow' of these couplings in the 
limit where the cut-off tends to zero, that is, as the abso- 
lute temperature vanishes. In the perturbative approach 
to RG, the flow is described by a set of differential equa- 
tions, whose solutions we study in this section. 



Simple power-counting arguments show that 



J dqdu q'^\uj\ |(/)(g,a;)p is an irrelevant perturbation in 
the renormalization-group sense. This is true provided 
D{q — 0) is not singular, which is indeed the case (see 
Fig. |4] and Appendix [b|) . Indeed, this term alone leads 
to a momentum dependent broadening of the long-wave 
length phonon excitations of the gapless phase of the 
model in Eq. 21). Therefore, in order to study the low- 



temperature properties of the model, it is justified to drop 
Sjy, and therefore we shall next focus our attention on 



the second term in Eq. (35) and consider the effective 
model described hy S = Sb + S^, where Sb is given by 
Eq. ^ and S"^ given by Eq. (|37|. In the half-filled case, 
we also have to take into account given by Eq. ([38l. 



The resulting action contains only marginal and (poten- 
tially) relevant perturbations in the RG sense, which we 
shall analyze in this section. In what follows, we shall 
consider the cases of integer and half-integer lattice fill- 
ing separately. The details of the perturbative derivation 
of the RG equations are given in the Appendix [Cl 



A. Integer Lattice filling 



To 0{gbD,9u) the flow equations in this case read: 



dgu 
dl 
dgtD 
d£ 
dK 

111 
dv 

di 



(2 - K)g^, 
il-2K)gbD, 
-igl + 2TTgtD)K\ 
-2'n:gi,DKv. 



(42) 
(43) 
(44) 
(45) 



We neglect terms of 0{gl^) or higher because gbfoiO) oc 
g%p, that is, gbo is already second order in the Bose- 
Fermi coupling, which is assumed to be small. For gto = 
0, the equations reduce to those of a pure ID Boson sys- 
tem in a commensurate potential flrst obtained by Hal- 
dane |48] (see also [44, 47j); for g^ — 0, the equations re- 
duce to those derived in Ref . [321 which describe the quan- 
tum phase transition between a Tomonaga-Luttinger liq- 
uid and a dissipative insulator (DI). 

The above equations show that near the SF to MI 
quantum critical point (corresponding to K* = 2, gu — 0, 
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gbD = 0) the dissipative interaction is a highly irrelevant 
operator because 1 — 2K w —3. Thus, the most im- 
portant effect of the Fermi component of the mixture is 
to introduce a rcnormalization of the periodic potential 
and the screening of the interactions, which leads to the 
rcnormalization of the Luttinger parameter K and the 
sound velocity v given by Eq. ( 33 ) . 



From the analysis of the RG equations, which implies 
that the dissipation is an irrelevant operator in the RG 
sense, we conclude that dissipative effects are weak in 
the MI phase where gu grows as the energy cut-off ™^ 

the absolute temperature) decreases. Thus, the dis- 
sipative term can be treated using perturbation theory, 
and leads to a small (when compared to the excitation 
energy) broadening of the phonon excitations in the su- 
perfiuid TLL phase. As for the excitations of the MI 
phase, which corresponds to a 'particle' (i.e. excess by one 
bosons) or a 'hole' (i.e. absence of bosons) propagating 
against the Mott-insulating background, the dissipative 
part of the interaction with the Fermi gas similarly intro- 
duces damping on their motion, which translates into the 
broadening of the excitation energy dispersion. Such en- 
hancement of the excitation broadening can be measured 
by lattice modulation spectroscopy [371 SSI ESI [M] • 



B. Half-Integer Lattice filling 

In this case, and given that the initial conditions are 
the same for the and 5*1, we note that they can be 
combined into a single term = S'[)[(j)]+S^[(j)], which 

can be written as: 

2ao 



(46) 

where goiO) = ^ [gboiO) + guoiO)]- The RG flow equa- 
tions for this system then read: 



dgu 
M 
dgp 

de 

dK 

~dl 
dv 



-{4gl + 2TTgbD)K^, 



-2T:gDKv. 



(47) 
(48) 
(49) 
(50) 



These RG equations describe the flow in the vicinity of a 
quantum critical point located at K* = g* = = 0. 
Integrating them numerically, we obtain the phase dia- 
gram depicted in Fig. [5) Thus, we find that, for a rela- 
tively weak boson-fermion coupling \gBF\/fJ-B ^ 10^^, 
the part of the phase diagram occupied by the SF 
Tomonaga-Luttinger liquid phase (TLL) shrinks consid- 
erably. The latter phase is identified by the RG flows 
for which both g^ and gjj — )■ as the phonon cut-off 
haoe^ is reduced to zero (i.e. for i — )■ -l-oo), that is, as the 
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FIG. 5: Phase diagram for the Tomonaga-Luttinger liquid 
(TLL) to Charge Density Wave (CDW) transition in the pres- 
ence of a Fermi gas for go — 2.5 x 10"*, which corresponds 
to IgBpl/fJ-B ~ 10~^ {fj,B being the chemical potential of the 
bosons). K is the Luttinger parameter of the bosons in the 
mixture (cf. Eq. 33 1 and gu oc Ub\\ + O^gBp), where Ub\\ is 
the external periodic potential. The shaded area is the TLL 
phase. The diagonal dashed line represents the TLL-CDW 
phase boundary in the absence of fermions. The curves in the 
diagram represent RG flows for the K and gu couplings for 
a set of initial conditions lying on the quarter circle on the 
right. The flow proceeds from right to left as as K always 
decreases according to Eq. |49[ 



absolute temperature is decreased. On the other hand, 
the CDW phase is identifled with those flows for which 
gu ^ 1 at a certain value of f . However, it is also worth 
noticing that we have observed numerically (see Fig. [6]) 
that, especially close to the phase boundary (red curve 
in Fig. [5]), gu{^*)/gD{(-*) ^ 1, even if 5„ becomes of or- 
der one first in all cases studied. This means that, even 
if the low-energy physics of this phase is dominated by 
the potential term oc g„, the dissipative effects are by no 
means negligible. It is interesting that this happens inde- 
pendently of how small the bare (?u(0) is, and even in the 
limit gjj(O) — )• 0+. This is because, ultimately, the RG 



fiow of gu {t) is controlled by the first term in Eq. ( 47 1 , 



which leads to a much faster growth, although for small 
f/„(0), the initial flow may be controlled by the second 
term in Eq. ( [iT] ). 

The RG flow equations indicate that the quantum 
phase transition occurs at K = 1/2, where the dissipation 
and periodic potential simultaneously become relevant, 
and the system is driven from a superfluid to a CDW 
Mott-insulating states. To study the interplay between 
the dissipation and interaction around the critical point, 
we adopt a variational self-consistent harmonic approxi- 
mation (SCHA) by choosing a trial effective action of the 
from: 



SM = 



dqdu 



i2n) 



(51) 



where we have deflned the Green's function Gv{q,uj) = 

-1 

with the dimension- 
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FIG. 6: Runaway renormalization-group (RG) flow of the 
couplings Qu, representing the periodic potential, and go, 
representing the efi'ect of the fermion-induced dissipation for 
K < ^. We found that even for relatively small initial poten- 
tial gu{0) the RG flow of gu{£) eventually overcomes the flow 
of goit) and becomes gu{^) ~ 1 first. This means that the 
system localizes and becomes a Mott insulator. However, as 
this plot illustrates, the effect of go , i.e. the renormalized 
dissipative coupling, is not negligible. 



FIG. 7: Mott gap A in the presence and absence of fermion- 
induced dissipation as obtained from the self-consistent har- 
monic approximation (SCHA, see Sec. IIIB for details). It 
can be seen that the dissipation greatly enhances the Mott 
gap by suppressing the quantum fluctuations of the bosons in 
the CDW Mott insulating state. Note that the SCHA erro- 
neous yields a discontinuous phase transition at the critical 



point K* = 
mation |T7| . 



This is a well known artifact of this approxi- 



less self-consistent parameters rj and A that can be deter- 
mined by the minimization of the variational free-energy. 
A variational estimate F^ar of the true free-energy F can 
be obtained from Feynman's variational principle [47) : 



IV. COMMENSURATE - INCOMMENSURATE 
TRANSITION IN THE PRESENCE OF 
DISSIPATION 



F < F 



(52) 



A. Integer filling 



Therefore, optimizing 6F^ar[Gv\/ SG^ = 0,) the param- 
eters 77 and A are found by solving the self-consistent 
equation above (Eq. (D2)), so that (see appendix [P] for 
further details): 



A 



8g« 



(2^)2 



a (77, A, X), 



7-lKTr+2VKTrA 



(53) 
(54) 

2K 



where we have introduced a(ry, A) 

The numerical solution of these equation for the gap A is 
shown in Fig. [7) It can be seen that the gap is enhanced 
for K < ^. This expected is because quantum dissipa- 
tion is akin to classical friction, which hinders the mo- 
tion of the particles and thus helps to stabilize the CDW 
Mott-insulating state. Note, however, that the SCHA 
erroneously yields a discontinuous transition for K = ^. 
This is a well known artifact of this approximation |47j . 



In this case, as for the TLL to Ml transition, the effect 
of dissipation is rather weak. A way of understanding 
this is to stop the RG flow when (?„(^) ^ 1 and consider 
the sine-Gordon model at the Luther-Emery point where 
it maps to a ID relativistic model of massive (Dirac) 
fermions lU 27]. Diagonalization of this model yields 
two bands separated by a gap: a filled 'valence' band 
and an empty 'conduction' band HH H?]- Tuning the 
chemical for the bosons amounts to introducing particles 
in the conduction band or holes in the valence band [JT] . 
For small particle (hole) density, the system can be de- 
scribed as a Tonks-Girardeau gas |5H] characterized by 
Luttinger parameter K 1. The dissipation being an ir- 
relevant for K > K* = i , its effect on such a dilute liquid 
of particles (holes) is negligible as far as the ground state 
properties are concerned (although it will lead to a small 
linewidth of the excitations, which is due to collisions 
between the bosons and the fermions). Thus, in partic- 
ular, the exponents characterizing the commensurate to 
incommensurate (C-IC) transition are thus expected to 
remain unchanged and, therefore, the density of particles 
(or holes) [HI will grow as -^/x — where fic ^, 
where A is the Ml gap. 
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B. Half- integer filling 

For half-integer filling the situation is very different, as 
it was already pointed out in our discussion of the previ- 
ous section. We can realize this by considering again the 
case where we take go infinitesimally small but ^ 1. 
Applying the same reasoning used in the previous section, 
the sine- Gordon model S}j[4i\ -\- Su[4>\ in this case maps 
to a system of Dirac fermions describing the (fraction- 
ally charged) soliton and anti-soliton excitations of the 
CDW state (configurations of the form 10101011010101, 
for the solitons, and 101010010101, for the anti-solitons). 
A dilute gas of such excitations can be described as a 
Luttinger gas with a parameter K :^ j < K* = Thus, 
the dissipative term Sd[4)\ from Eq. (46) is a strongly 



relevant perturbation, which, as discussed in Ref. [39] . 
leads to the localization of the system in a new phase, 
which we term dissipative insulator (DI). In this phase, 
the boson density, {pb{x)) exhibits long-range order [39] 
with a characteristic wave number equal to 47r/5o- 

However, it is worth mentioning that, as Fig.[6]demon- 
strates, assuming that is infinitesimal when go ^ ^ 
is not representative of the the RG flow described in the 
previous section. Indeed, we numerically found that even 
in the case 5„(0) 0, goi^*) < gu{t) ~ 1 (see Fig. [6]) 
in other words, the dissipation, although diverging less 
strongly than the periodic potential, is not a small per- 
turbation on the CDW state. Thus, we expect that the 
dissipative term needs to be treated on equal footing with 
the potential term oc g^- The universality class of the 
commensurate to incommensurate transition is therefore 
expected to be different from the case of integer filling. 



CONCLUSIONS 



dissipation is much more severe when the bosons are close 
to a superfluid to CDW Mott-insulator transition, which 
happens at half-integer filling. In this case, the dissipa- 
tive effects strongly hinder the motion of the bosons and 
help stabilizing the CDW phase (cf. Fig. [5| as well as 
enhancing the CDW gap (cf. Fig. [7|. This effect leads to 
a dramatic suppression of the superfiuid phase relatative 
to the pure boson case, which can observed as a reduction 
of the potential depth required for the bosons to local- 
ize in the CDW phase. The enhancement of the gap on 
the CDW side of the transition can be also probed using 
lattice modulation spectroscopy. 

We have also studied the commensurate- 
incommensurate transition and argued that in the 
case of integer lattice filling, the fermion-induced dissi- 
pation is an irrelevant perturbation and therefore, the 
universality class should not be altered. However, in the 
case of half-integer filling the dissipation is relevant (but 
less than the external potential) and therefore we expect 
the universality class will be modified. This subject 
requires further study, but it will not be pursued here. 
The conclusions of this work are summarized in the 
schematic phase diagram of Fig. [8) 
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In conclusion, we have studied a model for a mixed 
dimensional Bose-Fermi mixture in an optical lattice, 
where the bosons are confined to one dimension whereas 
the fermions are free to hop in three dimensions (albeit 
with renormalized dispersion) . We have argued that this 
system is a realization of a ID interacting Bose gas cou- 
pled to a dissipative bath of the Ohmic type. In addition, 
the fermions also screen the boson-boson interactions. 
For integer filling of the boson lattice, we have found that 
the dominant effect of the fermions on the bosons is the 
screening of their interactions, as it was also observed in 
mean-field studies of 3D dimensional optical lattices [25] . 
Thus, provided the so-called self-trapping effect can be 
subtracted or compensated, the screening of the boson 
interactions leads to an enhancement of the superfluid 
properties as the bosons become polarons with reduced 
effective interactions. In this case, dissipation effects only 
contribute to an increase in the linewidth of the excita- 
tions in both the superfluid and Mott-insulating phases, 
which could be detected by means of lattice modulation 
spectroscopy [5711151115] . 

On the other hand, the effect of the fermion-induced 




K=1/4 

FIG. 8; Schematic phase diagram of a system of heavy 
bosons confined to ID and and coupled to a dissipative bath 
of Ught fermions moving in 3D. 
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Appendix A: Relating x^{x,ijj) to xf{x,t) 



particle orbitals of the fermions are given by 



In this appendix we will derive the identity that we 
used in the main text to relate the retarded density cor- 
relation function to its imaginary time version at zero 
temperature. We shall first recall that the retarded cor- 
relation function is defined as: 



Xf (x,<) = --m {[6pF{x,t)5pF{QM 



(Al) 



e^Hpt/h^p^(^^y-iH^t/h^ and 5pF{r) = 



where 5pF{x,t) 
5pF{x,t) 

Pf{^) — However, the imaginary time correlation 

is defined as: 

Xf{x,t) - -^{SpF{x,T)SpFiO,0))F, (A2) 

where SpFix,T) = e" 5 pF{x)e-"''^'''. By taking the 
Fourier transform of the spectral representation of (Al) 
and comparing it to the spectral representation of 



Mj/2 

XF{x,iUn)= dTXF{x,T)e" 
J -hp/2 

we arrive at the following relation: 
du! Im x;P(a;, u) 



(A3) 



XF{x,iuJn) 



+°° dw 



(A4) 



ImxF(a;,w) ImxpiXyOj) 



UJ + iljJn 



(A5) 



where, in the deriving the last expression, we have used 
that Imx_F(a;, — w) = — Im xf{x, lu). Hence, introducing 
the last expression in (A3), taking /3 — )• +00, and per- 



forming the integral over uj^ with the help of Jordan's 
lemma, we arrive at the desired result: 



Xf{.x,t) = 



+00 



duj 



^1^1 Im XF^x^io). 



(A6) 



Appendix B: Fermion bath response function 

Let us consider the Fourier transform of the den- 
sity response of the Fermi gas at zero tempera- 
ture, which, as we neglect the interactions induced 
by the bosons on the fermions, is just the Lind- 
hard function. Recalling that the Matsubara ver- 
sion of the latter is defined as XF(a;, rj^, r^, t) — 
-\{5pF{x,Vi^,T)5pF{0,T'^,{)))F, where 5pf{x,y±_) = 
PF{x,r±) - p°{x,rj_), being pFix,r±) = 
Efc,fe',k^,k'^ '/'fc,k^(2;'i"_L)(/3fcMci(a;,r_L) fl^^fk'M'^ the 
density operator and p'^{x,r±) = {pFix,r±))F, the 
equilibrium density. We shall assume that the single 



ipk,\^^{x,r±) = ipk{x)ipw±{'r±) 
1 



R 



(Bl) 



where L is the (normalization) length in ID and M is the 
number of lattice sites labelled by R = (n, m)bo (bg is the 
lattice parameter), and WQ{r±) is Wannier orbital for the 
fermions. In the above expression we have assumed that 
the strength of the longitudinal potential in ID is weak 
so that the Bloch orbitals ipk{x) ~ Thus, we arrive 
at the following expression: 

XF(9,ri,rl,..)= J dre^(---«^) XF{x,r^,r'^,T) 



X Ai,k'_L(r_L,i"l), 



(B2) 



where the function ^k^.k'^ (r_L, r'_L) — 
V'k_L(''-L)V'k^('''-L)'^k'_L(r±)'/3*,^(r'^). The single- 

particle dispersion of the fermions is 

e(A:,k^) = e,i (fc) + e(k_L) 

2tj^ (cos kybo + cos k^bo) (B3) 



2m*, 



where we have assumed that the longitudinal dispersion 
is approximated by a quadratic dispersion characterized 
by an effective mass m*p niF and transverse motion is 
described by a tight-biding dispersion characterized by a 
transverse hopping t±. 

Indeed, the response function in which we are inter- 
ested is not the Lindhard function, |B2[ but the following 
integral of it: 

Xi=^((7,w„) = J drj_dr'j_Foirj_,r'j_)xiq,rj_,r'j_,ujn), 

(B4) 

where i^o(rj_,rj_) — \wo{r±)wo{r±)\'^ , where wo{r) are 
the Wannier orbitals for the bosons in the lowest Bloch 
band. Thus, in order to compute ( |B4[ ), we need to con- 
sider the following integral: 

J drj_dr'j_ Fo(r_L,r^) Ak^.k'^ (r_L, r'_L) 
^ ^.(k..R-kV.R') 

R,R' '' 



K(r^-R)]*<(r^-R') 



1 

M 



drj_ |wo(r_L)|^ U^(r_L)| 



A 

M2' 



(B5) 

(B6) 
(B7) 
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where we have approximated wo(r_L) 



-/(27r£|^) and <(r^) 

/ [2TTppj^) and assumed that <C ^f±j 

so that we can neglect overlap between the Wannier 
orbitals for R 7^ R'. In the above expression, 



j dri|u;o(r±)n<(ri)| 
Hence, 



1 



(B8) 



Xf(9,«w„) 



^'^^^fck k' «^"-e(fc + '7,k^) + e(A:,kj_) 

(B9) 



Next, we take the thermodynamic limit, transform the 
sums over fc, kj^,k'_L into integrals, and introduce the 
density of states of the 2D (square) lattice of tubes [52] , 



-K 



f-V 



2 



(BIO) 



where K{z) denotes the complete elliptic integral of the 
first kind and W = 4ij_. Thus, the retarded response 
function (obtained from XF{<1^ i^n) by means of analytic 
continuation where iujn — > uj^ = uj + iO"*") can be rewrit- 
ten as follows 

I-+W 

Xpil^^n) = A dede' p{e)p{e') 
J-w 
dk 



nk,6 — nk+q,e' 
2n huj+ +£-£'- eii (fc + g) + eii (fc) 



(Bll) 



r+W 

^ A dede' p{e)p{e') 
J-w 



-w 

dk 



1 



huj+ + e - e' - e\\{k + q) + e|| (k) 
1 

huj+ + £ - e' + e||(fc + (?) - e||(fc)_ 
A I dede' p{e)p{e') / —nk,e 



w 



2tt 



1 



+ e - e' - e|| (fc + q) + e|| (fc) 
•+)] 



w — > —a; 



(B12) 



At zero temperature ^ = 0(eF~e~e||(fc)), where e^? = 
^jLf{T = 0) is the Fermi energy (note that ep > —W 
otherwise there will be no fermions in the mixture). 

Let us first consider (minus) the imaginary part of 







/ de 


1 dkei^ 


l-W J 





2m% 



X P(e) 



p [ huj + e — 



p [hhj — e 



h^q^ 


h^kq 


2m*p 


m*p 




h'^kq\ 




m*p J 



*-2 2 

where we have set e||(fc + — ~ 



(B13) 

I h?kq 



The above expression can be used to obtain the (imag- 
inary part of the) response for arbitrary w. However, 
we are only interested in the regime of small w, for 
which we can expand p {huj ± E{k^ q, e)) = p{E{k, q, e))± 
p'{E{k,q,e))huj + --- (where E{k,q,e) = £ - + 
2kq)) and therefore, to lowest order in w, 

Im [— (g, w)] ^ Ahuj / de dk p{e) 



X p \e 



2m] 



-iq^ + 2kq) 



ep- e- 



fee 
2m*„ 



(B14) 



In order to perform the integration over fc, we define from 
the constraints imposed by the Heaviside step function 

in Eq. (|B14[), kpie) — y '^j^i^F ~ s) for e < ep, and 



note that 



dk d, 



4 



2m] 



-iq^ + 2kq) 



h^q J- 



+ kF{e) t2 , 

dkdkp{e-—{q' + 2kq)) 



h^q 



p e 



2m] 



-{q^ + 2kp{e)q) 



p{s-^iq^-2kpie)q)) 



Thus, the expression is simplified and only the integra- 
tion over e remains: 

M-x'^iq^oj)] c Ahoj[ " ^) y ^ deOiep ~ e)p[e) 

p{e-^y^2kpie)q)yp(e-^y-2kpie)q) 



(B15) 



This expression can be numerically evaluated (cf. Fig.|4|. 
However, for g — > 0, further analytical progress is possible 
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by noting that p{e)p'{e) = ^d[p{e)]'^ / de , and hence, 

Im [-XpiQ ^ 0,w)] ~ AhLj / dekpie) 

J-W "-^ 

(B16) 

From which, upon integration by parts, we obtain: 
+ [p(min{T^,ej.})]2fcf (min{M^,eF}) 



fi2 



= A/io;— 1^ 



kpie) 
de 



[P{e)] 



-w 



kpie) 



(B17) 



Hence, by direct numerical evaluation of the above ex- 
pression we see that it is not singular, which implies 
that the q ^ term (denoted Sjj in Eq. (36l), can 
be neglected. In general, using Eq. B15 to evaluate 
Im [— %|^(q,aj)] for finite g, we find it is also a nonsingu- 
lar function of q in the neighborhood of q = 2k p = 2ttp^ . 
The results of a numerical evaluation of the integrals in 
equations (B15) and (B17) are displayed in Fig.|4j 

Finally, the real part of the response function is given 
by: 



Re[xf(g,c^)] = A / dede'p{e)p{e') 



dk 
2^ 



nk,£P 



1 



/min[W.,e f] 
dep{e 
-w 

f'min[VF,eF] 

X / de'P 
-w 



hw + e - e' ~ t\\{k + q) + e|| (fc) 
dk / h^k"^ 



+ (cli — ^ — a;) 



27r 



2m*, 



huj + e-s' - en {k + q) + en (k) 



A 



dep{e' 



-w 

inin[H^,ej= 
-W 



de'P 



l-kF{6) 

Pie') 



E 



27r 

- + {uj ^ -uj) 



(B18) 



where we have introduced E — huj + e — (fc + g) + 
e||(A;). Furthermore, by using the well-known Kramers- 
Kronig relations that connect the real and imaginary part 
of any complex function which is analytic in the upper 
half plane: 



Re[G^(R = 0,£)] = -P 



de' Im[G^(R = 0,e') 



(B19) 



we have that: 

Pie') 



P / de' 



E-e' 



P / de'- 



-l/7r)Im[G^(R = 0,e')] 



E-e' 



= Re[G'\Ti^O,E)] 



(B20) 



Then, we can rewrite equation (B18I, so that: 

^min[H^,eF] 



Re[x^{q,oj)]=A 



dep{e) 



'^^(^^ dk 



2tt 



Re[G'"(R = 0,£;)] 



A 
2^ 



-W 



+kF{s) 



dk 



kp{£) 



de p[e) 



g [fujj ^ e - 



h q h^kq 
2m% m% 



g ( —hu + £ — 



ft^q"^ h?kq 
2171% m% 



(B21) 



where 5(e) — — Re[G-^(R = 0,e)] is the Hilbert trans- 
form of the density of states in a 2D square lattice mod- 
eled by a tight-binding approximation |52j : 



g{e) ^P de 



J Pie') 



2_ 



/C(^) for|e|>W^, 
lM/C(f) for|£|<I^, 
(B22) 



In particular, the static limit uj — Q reads: 

min{XV,ei^} 



Xsiq) = Re [xf(9,c^ = 0) 



dkq \ e — - — — 
-kAe) V 2m^ 



A 



fi?kq 



de p{e) 



(B23) 



Therefore, it is possible to perform the calculation of the 
previous expression. 

At low frequencies (hw <^ fis < £f) we shall approxi- 
mate the response function of the Fermi gas by the two 
first terms in the series about uj ~ 0, i.e. 



XF(g,w) ~ Xsiq) - iTrujD{q), 



(B24) 



where S{q) = Im [— XfI*?)] /i^'^) Finally, we make use of 
the spectral properties of relating the retarded response 
function to its analytical continuation to imaginary fre- 
quencies derived in the Appendix [A] 



xiq,^n) = - 



dijj Imx(g, 



In particular, the static limit w„ = corresponds to: 

dio Im x('?, w) 



(B25) 



Xsiq) = x(9,o) 



(B26) 
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Adding and subtracting the static part, 



= Xs{q) 

and recalling that 



did 



lmx{q,uj) 



1 



1 



TT Ul 



iuj„ 



lUJn — to 



(B27) 



27T 

Xs{q)6{T) - 

Xs{q)S{T) - 



rfwlnix(g,w) 






2tt uj 


l-oo 


iuJn — LU 



+ 00 



dw Im x{q, w) 



TT W 



(B28) 



Thus, upon performing the above integral over a;„ using 
Cauchy's theorem, the following expression is obtained: 



/"H-oo 7 

Xiq.r) = Xs{q)5{T) + / — e-Hlmxf(g,w) 



(B29) 



Hence, introducing Eq. (B24) in the expression above 

D{q) 



x{q,T) ~ Xsiq) S{t) - 



\r\ + ny 



(B30) 



The first term describes the short time behavior, which 
is dominated by screening, whereas the second term de- 
scribes the long time behavior, which is dominated by 
dissipation. 



Appendix C: RG analysis at half- filling 

In order to obtain the RG flow equations, we con- 
sider the functional integral representation of the par- 
tition function: 



Z(re) = J [#] 



-S14>] 



where 



(CI) 



(C2) 



^0 [(j)] being the Gaussian part of the action (the first term 
in Eq. 21). When writting (CI), we have made explicit 



the dependence of the partition function on the short- 
distance cut-off ao ~ vTc- Note, however, that (up to a 
multiplicative constant), the partition function is inde- 
pendent of the cut-off, and we will base our subsequent 
analysis on this fact. For a general perturbation S'int[0] 
we cannot compute the partition function exactly. Thus, 



we resort to a perturbative expansion of Z[{1 + 5£)tc] 
(where Si <^ 1) in powers of Sint- 



Z 



[(1 + Se)T,] = Zo[(l + S£)t,]{i - {SiM (C3) 



} 



(C4) 



To deal with this expansion it is convenient to define the 
normal ordered vertex operators: 



. g2pi0(x) . . 



1 



o2pi0(x) 



(C5) 



where x = (wt, x) the limit ao — t- is implicitly under- 
stood. Then, when inserted in an expectation value, we 
have the following operator product expansions (OPE): 



. g2ip0(r) .. g-2j:p0(r') 



1 



[1 + 2ip{v - r')V(/)(R) 



-2p2[(r-r')V0(R)]2 + ...] : (C6) 

. g2jp0(r) .. g2ip0(r') ._ ^2p^ K . g4-ip0(R) . ^ 

where r = (ur, x), R = (r — r')/2 V = ((1/w)9t, dx) and 
oq — VTc is a short-distance cut-off. Next, let us consider 
the partition function at the scale (l + (5^)ao, where 51 > Q 
and 51 < 1: 

Z[{1 + 5l)aa] = Zo[(l + 5l)a^][l - {Smt) + ^(^L) + ' ' ' } 

(C8) 



where 



Su[(l)] = ^^Z4K / dxdr : cos4(/)(r) :, 



(C9) 



S 



^0 
9D 



S(x - x') 

Dm = / ^i-^i-' 1^ ^,12 • cos20(r) : 



x:cos20(r'): (CIO) 
where we have normal ordered the vertex operators. 



1. First order terms 



Now, let us consider the first order term {Sint) 
{Su) + {Sd}: 



-{Su) 



9uil + 5l) 



7r[{l + 5l)ao 



2-iK 



dr(: cos40(r) :) (Cll) 



When compared with the same operator at the scale ao, 
we find that: 



9u{l + 5l) 



= gu{i) 



[{l + Sl)]^-^^ 
gu{l + 51) = guil)[l - {2 - 4K)5l], 



(C12) 
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which inmediately leads to the differential equation: 
dgu{l) 



dl 



(2 -4if )<?„(/) 



(C13) 



Next, we consider: 



goil + Sl) 



/ drdr'^fc^(: cos20(r) :: cos20(r') :> 

i|r-r'|>ao(l+5/) |r-r'|^ 

(C14) 

To bring this expression to a form which can be compared 
with the same expression at the cut-off scale qq, we first 
split the integral on r and r' as follows: 



drdr' ■ ■ ■ 



|r-r'|>ao(l+<5i) 

drdr' ■ ■ 

|r-r'|>ao 



drdr' ■ 



ao{l+Sl)>\r-r'\>ao 



(C15) 



Thus, from the first term int the right hand-side of the 
above equation, we have that: 

9d{1 + SI) 



Tr[{l + 6l)aoV-^^ 

,S{x — x') 



drdr' 



|r-r'|>ao 



r — r' 



-/(: cos2(/)(r) :: cos2(/)(r') :) 



(C16) 



Hence, following the same procedure as before: 



[(l + ,50]i-2^ 
dgoil) 



9d{1) 



dl 



= {l-2K)gD{l). 



(C17) 



Next, we take up the contribution from the second term 
in Eq. [CT5] 



9d{1 + 51) 



,1-2K 



, ,5{x- x') 
drdr -r tttt x 



ao(l+<5/)>|r-r'|>ao 



cos2(/)(r) :: cos20(r') :) 



goil + Sl) 



drdr' ~ ^'^ ^ 

24-2A' J^^^,^si)>\r-r'\>ao |r - r' |2+2^f 



1 - 2[(r - r')V<?!)(R) 



goil + Sl) 



2a, 



1~2K 



, ,S{x- x') 
drdr -r tttt x 



/ao(l+(50>|r-r'|>ao 

al''{: cos40(R) :) 
Introducing u = r — r' leads to 



(C18) 



goim 



du 



SiUx) 

1-4K I I,, 12 

an Jao(l+(5;)>|u|>ao 1^1 



SI 



1-4K 
On an 



(C19) 



Hence, the second term in Eq. (C18) yields: 

goim 



2-4K 



dr : cos 40(r) 



(C20) 



Therefore, the flow equation for (i.e. Eq. C13) must 
be modified to: 

guii + Sl) , 
- TigniljSl = guil) 



[il + SI) 

g^il + SI) - [1 + (2 - 4K)Sl]guil) + 7rgoil)Sl 
dguil) 



dl 



i2~ 4K)g^il)^7rgDil) 



(C21) 



Finally, it is necessary to consider the first term in 



Eq. (|C18|. To this end, wc need to consider the following 

Siux) 



integral with r — r' = u = (ur, ux) : 
goil)Sl 



1~2K 



dUT 



ao(l+(5;)>|u|>a(, 



\2+2K 



a„{l+Sl) 



,2K 



ao 

ao{l+5l) 



1 



,l-AK 



'(-) 

Van/ 



ao 



2K 



= 51 



(C22) 



51 I dxdT{idr<i 



(C23) 



Thus, a term of the following form is generated: 

^gojl) 

vil) 

Note that this term has the same form of the ^^"^^ idr'P)'^ 
operator in the Gaussian action, 5*0 [</>]. When re- 
exponentiated, we fined that, upon comparing with the 
same action at the cut-off scale as: 



1 



27rif(/ + 5l)vil + 51) 
Hence, 

1 



v{l) 



1 



2TTKil)vil) 

guil) 



(C24) 



Atv "^'-'' 51 

Kil + Sl)vil + Sl) 2TTKil)vil) ' vH) 
_d / J_\ _ ingo 
dl \Kv) ~ V 



(C25) 



Furthermore, the coefficient of / dxdridx't')'^ is not renor- 
malized: 



vil + 51) vil) dfv 



= 



(C26) 



Kil + 5l) Kil) dl\K. 

From these equations we can extract the RG flow equa- 
tions for k and v. 

1 dv 
Kdl 

1 dv 1 d / 1 \ 4:TTgo 
'K^'dl ^ vJAk. ^ 

1 du d / 1 

h v—\ — 

K dl dl\Kj 



V 

= At: gov 



(C27) 



(C28) 



Thus, adding Eq. ( |C27[ ) and Eq. ( |C28| : 

. d 



dl 



d / 1 



2T:go. (C29) 
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2. Second order terms 

After considering the first order contributions, we need 
to take up the second order: 



(Sim 



(C30) 



We do not consider terms of order gugo or gj^ because 
go ex. g'^p is already second order and ^ 1 is consid- 
ered smah. Thus, taking: 



21 (-5' 



gu{l + Sl) 



2\^^[(l + (5/)ao]2-4^ ; 
[ drdr'{: cos40(r) :: cos40(r') :) (C31) 

■/|r-r'|>ao(l+<5n 



Again, we spUt the integral as in Eq. (C15), which leads 
to: 



1/ 9u{l) 



I 

Jao(l 



ao(l+50>|r-r'|>ao 



drdr'{: ( 1 - 8[(r - r')V0(R) 



1 

„/|8A'-4 



dR(|VR0(R)|^ 



ao(l+5i)>|r-r'|>ao 

2 



(-) 



SK-4 



COS (f) + 



-glSl / dv{: (V0(r)) :) + 



(C32) 



Thus, we need to revise our previously derived equations 
for the renormalization of the Gaussian action parame- 
ters: 



1 



Hence, the RG flow equations for both K and v: 



d / 1 

JlVK 
d 

Jl 



Kv^ dl 
V = —2TTgoKv 



I dv 47r 

= —[9D 

V 



(C39) 
(C40) 



Finally, in the analysis of the second order contribu- 
tions, we need to consider the term: 



Oigugo) 



2gu{l + 5l) 



goil + Sl) 



Tr[aa{l + diy 



4-2K 



ao{l + Sl)] 



1-2K 



dridr2dr3 



S{xi - X2) 
|ri-r2|2 ' 



cos4(/)(ri) cos4(/)(r2) cos40(r3) :) 
(C41) 



where the star (*) under the integral means that: |ri — 
r2| > ao(l + SI), |ri - ra] > ao(l + 61), and |r2 - r^] > 
ao{l+Sl). Let us consider the contribution resulting from 
the OPE when ri r2 (or equivalently ri — > r^): 



: cos4(^(ri) :: cos4(?!)(r2) 



2|r- 



,./\4K 



: cos2(?!)(R) 



(C42) 



Hence, as the above factor of 2 is cancelled by the two 
possible contractions ri r2 and r2 ra: 



Oigugo) 



2gUl + Sl) 



goil + Sl) 



^[ao(l + '50]*"'^' [aoil + Sl)] 



1-2K 



dp- 



1 



ao(l+<50>|p|>oo IpI^^ 

X (: cos4(?!)(R) cos4(/>(r3) :' 
4g„(0.9D(0 



Six - X3) 



2TrKil + 5l)vil + Sl) vil) nvH) 
vil + Sl) 2vil)glil)Sl vil) 



2TrKil)vil) 
(C33) 



7T[aoil + Sl)]^-^^ 

drdr' ^1^-^41 {■■ cos2(^(r) :: cos2(/)(r') :) (C43) 
|r — r'l'^ 

Therefore, we obtain the following differential equation: 
goil + Sl) iguil)gDil)Sl goil) _^ 



2TTKil -f SI) 
Hence, 



2TTKil) 



dl V Kv I V V 
d /■ V 

JlKK 



Avg, 



1 d / 1 
vJl\K 

4lil<) 



1 dv At: , 



Kv^ dl 
1 dv _ 

kIi ~ 



[go 



gz 

n 



(C34) 

(C35) 
(C36) 

(C37) 
(C38) 



[aoil + 5l)Y-^^ 
gD_ 
dl 



,ao(l + <5/)]i-2if 
= (1 - 2K)gD + Agogu 



,1-2K 



(C44) 



dv 

dK 
"dl ' 
dgu 
dl 

dgp 
dl 



set of RG flow equations 


reads: 


—A-Kg^Kv 


(C45) 


---iAgl + 27rgD)K^ 


(C46) 


= (2 - AK)gu + Trgo 


(C47) 


= il - 2K)gD + Agogu 


(C48) 
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Appendix D: SCHA 



Therefore, requiring 5Fyar[Gv]/ 5Gv = yields: 



We have adopted a variational self-consistent harmonic 
approximation (SCHA) by choosing a trial effective ac- 
tion such as in Eq. (51 1 . To find the variational estimate 



of the free-energy we have to perform the averages of 
the effective action with respect to the trial effective ac- 



tion, by using Sq (Eq. pi}), SyjEq,. (|22j)) with p = 2 
for half-lattice filling, S'^Eq. (js?])) and (Eq. 
Thus, the variational free-energy Rvar that follows from 
Eq. (l52l) wiU be: 



+ T 
-T 



dqduj 



In Gv{q,uj) 



1 



2ttK ' V 
9u 



( ^Vsq ) Gv{q,uj) 



T 



ao 

9D 



dx I drdr' - 

|r-T'|>To 



T 
ao 

{Sv}v 



dx 



drdr' - 



-4/f^[l-cosc^(r-r')]G„(g,c.) 

|r-r'|2 

'4/f|^[l+cosa;(r-r')]G„(9,c^) 



\t-t'\>to 



r'|2 



(Dl) 



SFyar[Gv]/SGi, - [ r 

2 J (2tt)^ Gy{q,uj) 



T 



1 ,U} 



2ttK ' V 
(2^)2ao 



(— +w,g") 



(27r)2aoTo 



a(r7,A,AO + a2(?7,A,X) 



= (D2) 



where a{'q, IS., K) 



7-1 K IT +2 Kit A 



2K 



By keeping the 



r-independent terms in the integrals in Eq. ( |D1[ ) which 
yield the l eading contributions in lo to Gy{q,u;), leads to 
equations (53|54|. 
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